Excitation spectrum of Andreev billiards with a mixed phase space 
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We present a semiclassical theory for the excitation spectrum of a ballistic quantum dot weakly 
coupled to a superconductor, for the generic situation that the classical motion gives rise to a phase 
space containing islands of regularity in a chaotic sea. The density of low-energy excitations is 
determined by quantum energy scales that are related in a simple way to the morphology of the 
mixed phase space. An exact quantum mechanical computation for the annular billiard shows good 
agreement with the semiclassical predictions, in particular for the reduction of the excitation gap 
when the coupling to the regular regions is maximal. 

PACS numbers: 74.50.+r, 03.65.Sq, 05.45.+b, 74.80.Fp 



The spectral statistics of quantum systems is inti- 
mately related to the nature of the corresponding clas- 
sical dynamics jjj. Two celebrated examples are that 
chaoticity of the classical dynamics is reflected in the 
quantum realm by level repulsion while integrability 
causes level clustering pj. Recently, confined two- 
dimensional electron gases (quantum dots) coupled to a 
superconductor via a ballistic point contact have become 
a new arena for the study of quantum-classical correspon- 
dences [ pj-pl . Such systems are commonly called Andreev 
billiards M, because of the alternation of ballistic motion 
(as in a conventional billiard) with Andreev reflection ||] 
at the interface with the superconductor. The proxim- 
ity of the superconductor causes a depletion of excited 
states at low energies (proximity effect). It was found ||] 
that a chaotic Andreev billiard has an excitation gap of 
the order of the Thouless energy, while an integrable An- 
dreev billiard has no true gap but an approximately lin- 
early vanishing density of states. (The Thouless energy 
Et = gS/Ait is the product of the point contact conduc- 
tance g, in units of 2e 2 /h, and the mean level spacing S 
of the isolated billiard.) 

Both chaotic and integrable dynamics are atypical. 
The generic situation is a mixed phase space, with "is- 
lands" of regularity separated from chaotic "seas" by im- 
penetrable dynamical barriers. A generally applicable 
theory for the proximity effect in ballistic systems should 
address the case of a mixed phase space. In this paper 
we present such a theory. 

In a semiclassical approach we link the excitation spec- 
trum quantitatively and qualitatively to the morphology 
of non-communicating regions in phase space. Different 
regions exhibit greatly varying length scales, which also 
depend sensitively on the position of the point contact. 
Still, we find that a general relation exists (in terms of ef- 
fective Thouless energies) between these classical length 
scales and the corresponding quantum energy scales. The 
results for integrable and fully chaotic motion are recov- 
ered as special cases. For the mixed phase space our 
main finding is a reduction of the excitation gap below 
the value Et of fully chaotic systems. The reduction can 
be an order of magnitude, as we illustrate by a numerical 



calculation in the annular billiard M shown in Fig. 




FIG. 1. Andreev billiard consisting of a confined normal 
conducting region interfacing with a superconductor (shaded) 
over a distance W. The normal region is shaped like an annu- 
lar billiard, bounded by two excentric circles (outer radius R, 
inner radius r, distance of origins p). This figure represents 
the case R — 1, r — 0.35, p — 0.1, W — 0.8. A periodic tra- 
jectory is indicated, involving two Andreev reflections at the 
interface. For the Poincare map one monitors the collisions 
with the outer boundary (angle of incidence a and coordinate 
s along the boundary, with a — denoting normal incidence 
and s — denoting the point closest to the inner circle). 

We consider a two-dimensional ballistic region fa "bil- 
liard") of area A (mean level spacing <5 = 2ttH /mA) 
that is connected to a superconductor by an opening 
of width W (corresponding to a dimensionless conduc- 
tance g — 2W/Xf, where is the Fermi wavelength). 
Classical trajectories consist of straight lines inside the 
billiard, with specular reflections at the boundaries and 
retro-reflections (= Andreev reflections) at the interface 
with the superconductor. We assume that 5 <C Et <C A, 
where A is the excitation gap in the bulk superconductor. 
The first condition, 6 <C Et or W » Xp, is required for a 
semiclassical treatment. The second condition, Et <C A, 
ensures that the excitation spectrum becomes indepen- 
dent of the properties of the superconductor. 

The Andreev billiard has a discrete spectrum for e < 
A. (The excitation energy e > is measured with respect 
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to the Fermi energy. We count each spin-degenerate level 
once.) For e <C A the semiclassical expression for the 
density of states p(e) reads j| 
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with P(L) the distribution of path lengths between sub- 
sequent Andreev reflections. The distribution P(L) is 
normalised to unity and based on a measure ds d sin a 
for the initial conditions at the interface with the super- 
conductor (see Fig. |l|). The length scale Lt = hvp/2Ex 
(with vf the Fermi velocity) is determined by the ge- 
ometry of the billiard by Lt — irA/W and is therefore 
purely classical. Eq. ([!]) follows directly from the Bohr- 
Sommerfeld quantisation rule Q . The derivation of Lod- 
der and Nazarov || starts from the Eilenberger equation 
JTo| for the quasiclassical Green function and arrives at 
an expression for p{e) that is mathematically equivalent 
to Eq. (g). 

Return probabilities like P{L) and the related decay of 
classical correlations have been addressed in many stud- 
ies jllj . In a chaotic billiard, Lt is the mean path length 
and P{L) ex exp(— L/ Lt) is an exponential distribution. 
Eq. (nl) then gives the density of states 



scale Lt,i, the distribution Pi(L) decays exponentially 
for chaotic parts of phase space while algebraic decay is 
found for regular regions Q . In each case the partial 
density of states pi rises to a value 2 /Si on an energy 
scale i?T,ij but while pi has an excitation gap for the 
chaotic regions it rises linearly for the regular regions. 
Eq. (||) applies to those regions that are accessible for a 
given location of the interface. We call these "connected" 
regions. The other "disconnected" regions (usually some 
of the regular islands) do not feel the proximity of the 
superconductor and give a constant background contri- 
bution pi(e) = 2/5i in the semiclassical approximation. 

Two phase space measures Oi and Vi determine the 
mean length Lt.i = Vi/Oi between Andreev reflections, 
the effective level spacing 5i = (2ttK) 2 /mVi, and the effec- 
tive Thouless energy Et,i = Tivf /2LT,i = giSi/Air, where 
9i = Oi/Xf is the effective dimensionless conductance. 
The first is the area Oi that the region overlaps with the 
superconducting interface on the Poincare map (see Fig. 
H). It is a measure of the coupling strength of a region 
to the superconductor. The second is the volume Vi that 
the region fills out in the full phase space (r, (j>), where 
the coordinate r extends over the area of the billiard and 
4> G [0,27r) is the direction of momentum. 



2x coshx 
5 sinh 2 x 



ttEj 



(2) 



which drops from 2/5 (the factor of two arises because 
both electron and hole excitations contribute at posi- 
tive e) to exponentially small values as e drops below 
~ 0.5 Et- Eq. (0) is close to the exact quantum me- 
chanical result |||7 which has p = for e < 0.6 Et- For 
integrable motion P(L) decays algebraically cx L~ p with 
p close to 3. Eq. (jf|) then gives p(e) oc e p ~ 2 , hence an 
approximately linearly vanishing density of states. Nu- 
merical studies on the circular and rectangular billiard 
confirm the validity of the semiclassical approach ||[L2) . 

The semiclassical expression (Q) holds also for mixed 
dynamics. It allows to regard each non-communicating 
region in phase space as a distinct system, to be la- 
belled by an index i. It is helpful to rewrite Eq. ([!]) in 
terms of an effective level spacing Si and Thouless energy 
Eta = hvF/2LT,i for each of these regions. (This ap- 
proach extends the Berry-Robnik conjecture ]13[ to open 
systems.) We decompose p = J2i Pi m to partial densities 
of states pi, defined by 
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The distribution Pi(L) (still normalised to unity) 
now pertains to initial conditions (still with measure 
ds d sin a) on the interface with the superconductor that 
evolve into the i-th region of phase space. On the 
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FIG. 2. Poincare map of the annular billiard of Fig. [l]. Dy- 
namical barriers separating regions in phase space are shown 
as dashed lines. Chaotic trajectories are found in region 1. 
Region 2 is an island of regular motion around a short stable 
periodic orbit. Region 3 is integrable and consists of skipping 
orbits that never hit the inner circle. Initial conditions at the 
interface with the superconductor are uniformly distributed 
within the strip around s — tv, of area 2W . The hatched area 
is the overlap 0\ of this strip with region 1. 

The phase space that is explored from the point con- 
tact can also be parameterised by the variables s, sin a 
on the interface and a coordinate I along the trajectory. 
The identification of Lt.i — Vi/Oi as the mean path 
length in region i is a consequence of ds d sin a dl — dr d(f>. 
The mean length of all trajectories (L) = J dL LP(L) = 
J £2 i 'OiLT,i/2W — ^2iVi/2W can be used to determine 
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the total phase space volume V con = Yli^i = 2FF(L) 
that is connected to the interface with the superconduc- 
tor. Here the prime denotes restriction of the sum to con- 
nected regions, and we used the sum rule ^Z/Oj = 2W. 
The volume Vdi s of the disconnected regions (which de- 
termines the background contribution to p) follows from 
the sum rule £\ V< = V con + V&b = 2-kA. 

Since typically the smallest Er,i <C Et, the total den- 
sity of states p = ^\ pi has a reduced excitation gap. 
This is especially the case when one couples maximally 
to the regular regions. Then their contribution to p at 
small e (long path lengths) is minimal, and the gap is sub- 
stantially reduced due to long chaotic trajectories. The 
constant background and the linear increase from regu- 
lar regions dominates when the coupling is mainly to the 
chaotic parts of phase space. 

The preceding paragraph summarises the key finding 
of our work. We illustrate it now for the annular billiard 
of Fig. [P. T he Poincare map in Fig. || shows three main 
regions |15[, one with chaotic motion (1) and two with 
regular motion (2 and 3). The regular island 2 corre- 
sponds to orbits that bounce back and forth between the 
two circles where their distance is largest. It has a short 
stable periodic orbit in its centre. Region 3 is integrable 
and consists of skipping orbits (trajectories that do not 
hit the inner circle, so that their angular momentum sin a 
is conserved). 
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FIG. 3. Density of states of the annular billiard of Fig. |l| 
The solid curve is the semiclassical prediction computed from 
Eq. ([!]) . The histogram is obtained by an exact quantum me- 
chanical computation. The dashed curve is the semiclassical 
result ^ for completely chaotic dynamics. 

The regular regions couple maximally to the point con- 
tact when it is located at the short stable periodic orbit, 
as in Fig. [I] (location s = 7r). We have computed P(L) 
by following trajectories and obtained p(e) from Eq. (Q). 
The result is shown in Fig. [| (solid curve). We see an ex- 
citation gap which is about a factor of 4 smaller than in 
the fully chaotic case [Eq. (|^) , dashed curve] . The reduc- 
tion originates from long chaotic trajectories with mean 



length Lt,i ~ 4Lt, hence Et i ~ -Er/4. An exact quan- 
tum mechanical calculation flq| (histogram) confirms the 
low-e behaviour found semiclassically. The sharp peaks 
at higher e in the semiclassical prediction, which arise 
from families of regular trajectories of almost identical 
length, are not resolved in the histogram. This is not a 
surprise since numerically accessible Fermi wavelengths 
are still larger than the extension of the families. 
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FIG. 4. 
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Density of states of the billiard of Fig. | 
with two different locations of the interface [s — in(; 
s = 1 in (b)]. The semiclassical prediction from Eq. (Ill) (solid 
curves) is compared with an exact quantum mechanical com- 
putation (histograms). 

The regular island is disconnected from the supercon- 
ductor when the point contact is moved to the other end 
of the billiard (at s — 0, where the separation of the cir- 
cles is smallest). The gap in the chaotic partial density of 
states is reduced to a lesser degree than before, see Fig. 
||(a). Excitations localised in the regular island give a 
constant background contribution 2/8 2 = 2mV2 / (2-kTi) 2 . 
If the point contact is placed between these two extreme 
positions (at s — 1), the regular regions of phase space 
dominate the low-energy behaviour of p(e). Instead of an 
excitation gap we observe a smoothly and slowly increas- 
ing density of states, see Fig.|(b). Thehisto grams in Fig. 
^ fall systematically below the semiclassical prediction. 
We attribute this discrepancy to the constant background 
contribution in the semiclassical result, which should van- 
ish at small e because of quantum mechanical tunnelling 
through the dynamical barrier between regions 1 and 2. 
This source of error is absent in Fig. |^, because there all 
regions are directly coupled to the superconductor. 
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In summary, we have found that the superconductor 
proximity effect in ballistic systems depends sensitively 
on the morphology of the classical phase space. The ex- 
citation spectrum at low energies can be described in an 
intuitively appealing way by means of effective Thouless 
energies and level spacings for the regular and chaotic 
regions of phase space. If the coupling to the regular 
regions is maximal, the excitation spectrum exhibits an 
excitation gap that is much smaller than the gap of a 
fully chaotic system. Measurement of such a reduced 
gap would provide a unique insight into the effect of a 
mixed classical phase space on superconductivity. 
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